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Perturbative QCD results on tt production are presented in proton-proton, proton- 
nucleus and nucleus-nucleus collisions from CERN SPS up to RHIC energy. A 
Kj e t(s,pT,Q) factor obtained from jet production is applied to perform next- 
to-leading order calculations. Using the intrinsic transverse momentum (fcr) we 
determined transverse momentum spectra for pions in wide energy region. We 
have investigated nuclear multiscattering and the Cronin effect at RHIC energies. 



1. Introduction 

The calculation of high-pT particle production at SPS and RHIC energies 
requires solid theoretical background, based on perturbative QCD (pQCD). 
In these calculations, the choice of renormalization, factorization and frag- 
mentation scales are typically chosen in the range pt/3 < Q < 2pr- Earlier 
we have found 1 that the theoretical reproduction of the measured nuclear 
effects in pA and A A collisions (specifically the Cronin effect 2 ) strongly de- 
pends on the choice of the above scales. In the following we display new re- 
sults analysing preliminary RHIC data 3 on pion production in AuAu — > ir° 
at y/s = 130 AGeV and 200 AGeV. 



i 
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2. Perturbative QCD calculations with intrinsic kj> 

In pQCD-improved parton model (including intrinsic fcr) pion production 
was calculated in leading order (LO) 4,5 and next-to-lcading order (NLO) 6 . 
An intermediate solution to speed up the time consuming full NLO calcula- 
tions is to determine a Kj et (s,pT,j e t, Q) factor corresponds to higher order 
contributions at jet level in pp collisions 7 and apply it in the factorized 
pQCD equation. In this case the pion production in AA collision can be 
described as: 

K-^- = ^2 d 2 bd 2 rt A {r)t B (\b-r\) / d 2 fc T , a d 2 fc T , b dx a dx b dz c x 

" P abed J J 

X 9pA{k T .a) fa/A(Xa,Q 2 ) g P A{k T ,b) /b/A^b, 0?) X 

- ^ s d(s + t + u) . (1) 

dt 



ab—>cd 

K jet (s,p T ,jet,Q) — 



The partonic cross section, da/dt, is the LO Born term, which is multiplied 
by Kj et (s,pTjet,Q) m e Q- (!)• (s, i and u are parton-level Mandelstam 
variables.) The nuclear thickness function, = J dz p(b, z), is normal- 

ized as J d 2 btA(b) = A. The nuclear parton distribution functions (PDF) 
fa/A{ x ,Q 2 ) are based on the nucleon PDFs (f a / P (x a ,Q 2 )) and modified 
in the nuclear environment ("shadowing") 8 . We are using NLO MRST 9 
PDF set with Q = Prjet scale, where Prjet — Pt / z c- 

In a phenomenological approach, we introduced an extra 2-dimensional 
Gaussian transverse momentum distribution as the PDF's extension 5 : 

9 P A(k T ) = ,, 2 \ e - fe T/< fc T>P.4, where (k^) pA = {k 2 T ) pp + C -h pA {b). (2) 

We have found the Gaussian width ((k^pp) to be pt independent in the 
window 2 GeV < pr < 6 GeV. Nuclear multiscatterings yield an extra 
broadening in pA collision, which can be related to the number of nucleon- 
nucleon (NN) collisions in the medium and denoted by C-h p A(b). Here 
C means the average increase in the transverse momentum width in one 
NN collision and h p A{b) is a geometrical effectivity function. This latter 
quantity can be saturated 4,5 . The D^/ C {z c ^ Q' 2 ) is the fragmentation func- 
tion (FF), which gives the probability for parton c to fragment into 7r at 
momentum fraction z c and fragmentation scale Q' . In this paper we fix 
Q' = Pt/2 and use the KKP parametrization 10 . 
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3. pQCD result on pp and AuAu at RHIC energies 

In Figure 1 we display our results at y/s = 130 AGeV (left panels) and 200 
AGeV (right panels). In the upper panels the dashed curves indicate the 
pp results which can be compared to the preliminary PHENIX data 3 (full 
triangles) at \fs — 200 GeV. Full dots indicate the 7r° data from peripheral 
AuAu collisions, the calculated pion spectra (solid lines) approximately 
overlap with the data. We obtained =19 binary collisions for 60 — 

80% peripheral case and N u °® p — 10 for 70 — 80% peripheral case. A 
satisfactory agreement between the data and the calculations demonstrate 
the validity of the improved pQCD calculation^ 

Full squares indicate the data for central collisions. The number of bi- 
nary collisions are N M 3 ® C = 857 for 10% central AuAu case at yfs = 130 
AGeV, and N$° c = 940 at y/s~ = 200 AGeV. One can find a large differ- 
ence between central data and theoretical spectra (full lines), a suppression 
factor of <~ 4 — 5 appeared. 

Lower panels display the nuclear modification factor, defined as 

(Nun) <rf n A (b) dapp/dV 
This expression magnifies the enhancement of particle production in AuAu, 
relatively to the N^n upscaled pp collision on a linear scale. Applying 
shadowing 8 and multiscattering, the Cronin peak can be seen clearly in 
both peripheral (dotted lines) and central (solid lines) collisions. The posi- 
tion of Cronin peak depends on the {k^) pp value, and the height depends 
on the extra broadening 1 , C-h p A{b). One can see that the height of the peak 
in central case is larger than in peripheral one, which effect is caused by the 
number of nucleon-nucleon collisions inside a nucleus. We have fixed the 
PDF scale Q — prjet and FF scale Q 1 — pr/2 to obtain the peak's max- 
imum at pt ~ 3 GcV/c. The comparison of PHENIX data (full squares) 
with theoretical calculations display a difference of <~ 5. This suppression 
value is indicated on lower panels with a dashed lines. The origin of this 
suppression is jet quenching 11 . 



n AA ~ /AT __ \ _AA(U\~ 713— ■ y 6 > 



"Parameters: (k^)™ = 1.8 GeV 2 , = 1.5 GeV 2 , C sa t = 0.4 and v max = 3. 
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Figure 1. Kjet based NLO pQCD n° p^-spectra (upper panels) and the nuclear modifi- 
cation factor R\ A (lower panels ) compared to preliminary PHENIX AuAu and pp data 3 
at Vs = 130 AGeV (left) and 200 AGeV (right) RHIC energies. (See details in text.) 
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